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1. Introduction 

In the years following the discovery of Hawking radiation and the associated evapora- 
tion of black holes , there have been many efforts to either prove or refute the resulting 
implication that an initially pure state can collapse into a black hole and evaporate into a 
mixed state. The fact that such efforts have not proven successful is due to a combination 
of complications, including principally those of the backreaction of the Hawking radiation 
on the metric, and of the regions of large curvature (and hence strongly coupled quantum 
gravity effects) which are expected in gravitational collapse. 

Recently, Callan et. al. (CGHS) proposed a model which seemed to avoid some of 
these difficulties @. It consistts of gravity coupled to a dilaton and conformal matter in 
1 + 1 dimensions. For a single matter field it was found the backscatter (i.e., the Hawking 
radiation) occured in a region of strong coupling. By proliferating the number N of matter 
fields, it was believed that the essential physics was occur in a region of small coupling and 
hence be amenable to a systematic 1/N semiclassical expansion. 

These initial hopes were dashed 00 by the observation that the dilaton develops 
a singlurity at a finite value, dependent on N, precisely in the region where quantum 
fluctuations begin to become large. As a result, a number of groups have recently 

tried to explore, both numerically and analytically, the solutions of the large iV field 
equations. In particular, one is interested in the final "endpoint" of the Hawking radiation. 
Therefore, in ]|] || 0the fields were assumed to depend only on a "spatial" coordinate (of 
which there are a few natural choices). For example, in 0, a series of solutions with 
finite ADM mass and vanishing incoming and outgoing flux were found. Starting at weak 
coupling at spatial infinity, they were found to "bounce" back to weak coupling in the 
region of the singularity mentioned above. 

The static approximation used to derive these results is a significant simplification, 
but makes it difficult to consider the approach to the endpoint of the Hawking process. 
In the following, we will consider time-dependent (approximate) solutions to the CGHS 
equations. We will find solutions which still have finite ADM mass and vanishing flux, 
as well as regions with a time-dependent singular event horizon. In a later section, we 
will also discuss a series of perturbative, time dependant solutions which lie entirely in the 
Liouville region, followed by some concluding remarks. 
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2. The CGHS Model 

The CGHS model of dilaton gravity coupled to iV conformal matter fields in 1 + 1 
dimensions with coordinates a and r is defined by the action 

1 f 1 N 

S = — J d 2 aV^[e~ 2 HR + 4W) 2 + 4A 2 ) - - $>/,) 2 ], t 2 ' 1 ) 
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where g,(fi, and fi represent the metric, dilaton, and matter fields, respectively, and A is 
the cosmological constant. Integrating out the matter fields and going to conformal gauge, 
where 

9+- = ~\e 2p , (2.2) 
g++ = = 0, (2.3) 

(x± = r ± a), the resulting action is 

S =IJ d 2 a[e- 2 ^d + {24> - p)d.(2<p - p) - \ 2 e 2 ?) + - e^d+pd.p]. (2.4) 

The equations of motion for p and (p are 

T + _ = e- 2 *(2d+d-(j) - 40+00_0 - \ 2 e 2p ) - —d + d_p = 0, (2.5) 

-4<9 + <9_0 + Ad+<j)d_<i) + 2d+d-p + X 2 e 2p = 0. (2.6) 
Since the gauge has been fixed as in ( |2.3| ), there are two constraint conditions, namely, 

T±± = e- 2 *{Ad ± 4>d ± p - 281(f)) ~ ^(d ± pd ± - d\p + t±(a ± )) = 0, (2.7) 

where the functions t± are fixed by boundary conditions. 

The simplest and most important nontrivial solution of ( |2.5| )and ( |2.6[ )is the linear 
dilaton vacuum 

p = 0, cj> = ~(a+-a-). (2.8) 



2 

This vacuum has a singulrity at 



1 N 

--In—, (2.9) 

2 12' v ! 



as seen by calculating the sign of the kinetic operator in ( |2.4| ). As in previous papers, we 
will call the region of (j) < <J) cr the dilaton region, and (fi > (fi cr the Liouville, or strong 
coupling region. 



2 



3. Finite Mass Solutions 

Solutions to ( |2.5| ) and ( |2.6| ) with finite ADM mass were first found in 0by assuming 
that both (f> and p are time independant. In that case, (|2J^) and (|2.6f) become 

T+- = e- 2 \-\^" + 0'2 - X 2 e 2p ) + ^p" = 0, (3.1) 

4>" - <P' 2 - \p" + X 2 e 2p , (3.2) 

where the primes denote d/da. Linearizing about the linear dilaton vacuum solution (|2.8|) , 
for vanishing incoming and outgoing flux t±, asymptotically the resulting equations can 
be expressed as 

2\5<j)' + 2X 2 Sp - X5p' = 0, (3.3) 
5<P" = 8p"{l - ^e" 2 ^). (3.4) 
The asymptotic form of the solutions of these equations is 

^ = -^e- 2A ^ + ..., (3.5) 

5p = -^e-^ + ..., (3.6) 
where the parameter M is the ADM mass, given by evaluating 

M = 2e 2Xa (X5p + 5(P'). (3.7) 

at spatial infinity. 

Before going beyond the static case, it should be noted that one can expand 5<fi and 
SX in powers of e = e _2Acr , 

oo oo 

50 = ^a n e", 5p = Y,bne n , (3.8) 

n=l n=l 

with ai = bx = — M-. Substituing into the full linearized equations, one finds the relations 



1 - n 2 )a n+1 - b n+1 + 7^n 2 b n = 0, (3.9) 
(n + l)b n = 2na ni (3.10) 
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from which one easily finds 



jV n 3 

a n +i = — 3 2 ra n . (3.11) 

12 n d + — ?i — 1 



For large n, this suggests that we must have < 1 for the series to converge. For example, 



12 

for j| = 1, the resulting series for 8(f) is roughly 

M 1 



2Ae 



2Ac 



(3.12) 



thus implying that our linearizing approximation is breaking down for small a. It is perhaps 



of interest that the requirement ^ < 1 implies ( |2.9[) that (/> cr > so that the effective 
critical coupling constant e 2 ^ cr > 1. 

Let us now proceed beyond the static limit, but continue to require a finite ADM 
mass. Including time derivatives, in the a, r coordinate system the linearized equations 
read 

e 2p {~\H" ~ 2A5</>' - 2\ 2 5p + ±54>) + ^(5p" - 5p) = 0, (3.13) 
8(f)" - 6(j> + 2\8(f)' + -8p - -8p" + 2X 2 Sp = 0. (3.14) 



From (|3.7|) , we see that finite ADM mass requires both Sp and 8<f> vary asymptotically 



as e 2AcJ , as in (|3.5|) and (p.6|). If we express the perturbations about the linear dilaton 
vacuum as 

8(f) = x(r)e- 2Xa , (3.15) 

5p = y(r)e- 2X ^ (3.16) 
then to leading order ( |3.13j ) and ( |3.14| ) become, respectively, 

2A 2 x - 2X 2 y + -x = 0, (3.17) 

\y-x = Q. (3.18) 

It is a simple matter now to assume that x and y both vary as e WT and solve for u and the 
relative amplitudes. Of course, one solution is just 

M 

x = y=— + ar (3.19) 



as in (|3.5| ) (where a = 0). The other solution is easily seen to be 



x(t) = ^y(r) = ae~ 2XT + be 2Xr . (3.20) 

Li 



Substituting ( |3.20| )into (|3.7| ), we seen that the time dependancy of 5(f) and of Sp cancel, 



and the ADM mass is constant, even thought the metric and the dilaton are certainly not. 
Presumably, we should set the coefficient b in (|3.20 ) to zero, so that the solution is well 



behaved as r — > oo, as should the coefficient of the linear term in the uj = solution. 

The behavior of these solutions can be understood in much the same manner as in the 
static case ||. Let us concentrate on the M = solution, as it has been suggested that 
it represents the true quantum vacuum of the theory In any case, for r sufficiently 
negative, the time dependant terms dominate over the static terms. As one integrates the 
equation of motion in from spatial infinity, the solution may approach the singularity at 
(p cr (in the static case, this approach was guaranteed). In this region, we can essentially 
set p = 0, and cf> = <p cr + cp. The resulting equation of motion is 

-<p(<p"-<p) = ~(<p' 2 -<p 2 -\ 2 ). (3.21) 



If we continue to assume that (p = — 2A<^>, then (|3.21|) can be integrated, yielding 



*' 2 - - ~ ^AV = A 2 , (3.22) 

ip 3 

where A is an integration constant. As long as A ^ 0, this is the equation for a particle 
in a potential with an infinite barrier at the origin, so (p will bounce back to the weak 
coupling regime. 

We can also discuss the behavior of the solutions for any region where p — > — oo, in 
particular as a — ► — oo, assuming that ae~ 2Xr < as was discussed in the static case in 
0, by dropping terms proportional to e 2p which become irrelevant for p — ► — oo. For in 
that case we have 

e~ 2 ^ = -^-p + a+a+ + a-a-+b 1 (3.23) 



e- 2 * " ^ " ^H\j - g + e"*] = /(*_) + g(a + ), (3.24) 

where a± and b are constants, and / and g are arbitrary functions of their arguments (in 
the static case 0, one has f + g = —aa + c), the only priviso being that / must be smooth 



(i.e., /(cr_) is the integral of a completely arbitrary function). Concentrating on a region 
where / + g — > oo, we have 

e" 2 ^~ f + 9 + ^Hf + g), (3.25) 
12 1 

Using the formula for the curvature, 

R = 8e- 2p d + d_p, (3.27) 

we have 

R ~ d +9 d ~f e %{f+g-a + <T + -a-*--b)^ ( 3 28 ) 

(where we have redefined the constants a± and b). Taking, for example, g(cr+) ~ (cr+ — 
(T° ) _a , a > 0, we see that o+ is a singular event horizon. Since a+ = r + a, the location 
of the horizon is not constant in time r. Furthermore, the fact that bp grows more rapidly 
than S(f>, as seen in ( |3.20| ), suggests that such regions might be of greater importance in 
understanding the full evolution of the system, particularly for the M = solution, which 
has been proposed to be the true vacuum of the theory. In fact, in the original, unperturbed 
field equation ( |2.5| ) , we see that if <$>' = <j> = —2Xcf>, then p is forced to approach — oo, 
unless e 2 ^ ~ 24/N. Of course, at this point, depending on N, we may no longer be in the 
weak coupling regime which we have been discussing, but rather in the strong coupling, 
or Liouville region, which we consider below. 

Of course, for large r, the time dependant terms are small, and the solution behaves 
as in the static case, where <p penetrates closer and closer to <f> cr before bouncing back to 
weak coupling @]. But for r sufficiently large and negative, we are effectively dealing with 
the M = solution, in which p will tend to grow faster than <p and singular event horizons 
should appear. It is questionable whether or not this is a reasonable condition for the 
true vacuum of the theory. Actually, it seems more reasonable that the final state of the 
system, in response to some incoming matter, would have a potentially complicated causal 
structure. Of course, our solutions are nonsingular at 4> cr whereas the incoming matter is 
singular there, so the interpretation of these solutions remains unclear. 

To complement these solutions, we should in prinicple search for time dependant 
solutions with regular horizons, as was done in 0, H, and H], generally by using the 
"spatial" variable s = x+X- and then imposing continuity conditions at the horizon at 
s = 0. Including time dependant terms, of course, will affect the location of the horizon in 
general, and we have not yet made a determined effort to analyze the range of possibilities. 
Work on this problem is in progress. 
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4. The Liouville Region 

As argued in |7]], ||, solutions which lie entirely in the Liouville region contain im- 
portant information concerning the behavior of extremal four-dimensional dilaton black 
holes. Secondly, it might be possible that a configuration in the Liouville region might 
evolve into the weak coupling region, even if the reverse is impossible. 

To analyze this region, we introduce the new dependant variable 

iP = e~t (4.1) 

in terms of which the action is just 

S = - [ d 2 a(4d+ijd-ij + Atjjd+tpd-p - \ 2 ^ 2 e 2p + —d+pd-p). (4.2) 

TV J 12 

The resulting field equations (which can just as easily be derived from the original field 
equations upon substituting (|4.1f ) ) are 



T + _ = -2d+ipd-ip - 2^d + d-tP - X 2 tjj 2 e 2p - ^zd+d-p = 0, (4.3) 

4d+d_i; + 2i;d + d-p + \ 2 ^e 2p = 0. (4.4) 
The simplest solution to these equations is the trivial solution 

^ = 0,p = 0. (4.5) 

If we now perturb these equations about Q4.5Q , we see that every term in ( [L3|) is quadratic 
except the last term, so we just have 

5p = f + (a + ) + f-(a-), (4.6) 

where f± are arbitraty functions. Similarly, the linearization of ( |4.4| ) yields simply the 
Klein Gordon equation 

d+d-5tjj + —5iP = (4.7) 

for a particle with m 2 = A 2 /4. 

Another solution of @ (O) is § 



p = -ln(v^Ao-), (4.9) 

which is an example of anti-deSitter space, as the curvature turns out to be R = —4 A 2 . 
Linearizing again, we find 



IN N N 1 IN Sib 

-2\ —d-d + 5ip — —d-d + 5tp — — —^5p — \/ — — 5- = 0, (4.10) 
w 24 + ^ 12 + ^ 24 a 2 H V 24 <r 2 



and 

[W N N 1 

4 V 24 9 - a+ ^ + 12 a " a+ ^ + 24 ^ 6p = °" (4 ' n) 

Adding the equations, we have 

28-d+Si; - 5 4t = 0, (4.12) 

which is just the equation for a particle in a 1/r 2 potential. For example, going to the 
static limit, we have 

&if>" = ~5if>, (4.13) 
a z 

with solutions 

= a^ 1 +a 2 a 132 , (4.14) 

where the are constants and the fa are the solutions of the quadratic equation x 2 — x+2 = 
0. Since the fa are therefore complex, whereas ip should be real, it would seem that this is 
an inappropriate background for such a perturbative analysis. 



5. Discussion 

Spurred on in part by recent advances in string theory , we have witnessed a 
great increase in the number of toy models, particularly in low dimensions, made available 
for the study of phenomena such as Hawking radiation and the final state of black holes 
which involve fundamental issues surrounding quantum gravity. The CGHS model is an 
especially simple yet sufficiently rich example of such a model. Unfortunately, there remain 
significant barriers which interfere with our greater understanding of quantum gravity. Of 
the various groups who have studied the CGHS system, there are adherents of a variety 
of scenarios, including naked singularities [|J , macroscopic objects |l(| , the "bounce" 
scenario 101 , and so on. 



In this letter, we have tried to begin the program of going beyond the static limit 
applied earlier @ H @. We know that the classical no-hair theorems, which essentially 
say that a black hole is characterized by the quantum numbers of long range fields, such as 
mass, charge, angular momentum, cannot contain quantum mechanical information. What 
we have found is that specifying the mass of the black hole does not fully specify the metric 
or dilaton, even to leading order asymptotically There is active research underway on a 



variety of quantum-mechanical effects on black holes, see [lljfor example for a thorough 
discussion of quantum hair and Aharonov-Bohm type interactions of black holes. 

In the present case, in the original CGHS model (i.e., N = 1), the picture of the black 
hole was of an asymptotically flat plane connected via a throat-like horizon to a semi- 
infinite cylinder-like region. When matter impinges on this system, one might imagine, for 
example, that while the asymptotically flat region would eventually see a constant mass, 
the matter might be hurtling down the cylinder behind the event horizon in a complicated 
and possibly singular fashion. Even the horizon itself need not be fixed, though of course 
that would be measureable to an asymptotic observer. 

Another important factor which we have come across is the problem of the crossover 
between weak coupling and Liouville regions. In spite of the initial hopes, it appears that 
the important physics is occuring precisely in this region, where we cannot ignore futher 
quantum corrections. This region is small (of order A -1 ) in the large iV limit, so the model 
may yet be viable for questions regarding longer range phenomena. Furthermore, because 
of this great uncertainty, we cannot say for certain that propagation through the apparant 
singularity is in fact forbidden. Perhaps a further exploration of the appropriate boundary 
conditions or additional terms in the 1/N expansion will suggest a way out of our present 
dilemnas. 

Acknowledgements: The author would like to thank J. Lykken, S. Chaudhuri, H. Dyk- 
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